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Dark energy and dust matter phases from an exact /(i?) -cosmology model 
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We show that dust matter-dark energy combined phases can be achieved by the exact solution 
derived from a power law f{R) cosmological model. This example answers the query by which a 
dust-dominated decelerated phase, before dark-energy accelerated phase, is needed in order to form 
large scale structures. 
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^ Very recently, alternative theories of gravity are playing an interesting role to describe the today observed Universe. 



Although being the best fit to a wide range of data pj , the ACDM model is affected by strong theoretical shortcomings 



"^ \ [2] that have motivated the search for alternative models [3, 4]. 

QQ . Dark energy models mainly rely on the implicit assumption that Einstein's GR is the correct theory of gravity 

p\| ' indeed. Nevertheless, its validity on large astrophysical and cosmological scales has never been tested but only 

assumed [5| , and it is therefore conceivable that both cosmic speed up and missing matter are nothing else but signals 

of a breakdown of GR. In this sense, GR could fail in giving self-consistent pictures both at ultraviolet scales (early 

O ' universe) and at infrared scales (late universe). 

I Following this line of thinking, the "minimal" choice could be to take into account generic functions f{R) of the 

Ricci scalar R. The task for this extended theories should be to match the data under the "economic" requirement that 
no exotic dark ingredients have to be added, unless these are going to be found by means of fundamental experiments 
[a,[3|. This is the underlying philosophy of what are referred to as f{R)-gravity (see |J,l3)13 ^^-'^ references therein). 
Although higher order gravity theories have received much attention in cosmology, since they are naturally able to 
give rise to the accelerating expansion (both in the late and in the early universe 7]), it is possible to demonstrate 
that f{R) theories can also play a major role at astrophysical scales. In fact, modifying the gravity Lagrangian affects 
(^ • the gravitational potential in the low energy limit. Provided that the modified potential reduces to the Newtonian 
^H [ one on the Solar System scale, this implication could represent an intriguing opportunity rather than a shortcoming 
C^ ' for f{R) theories. In fact, a corrected gravitational potential could offer the possibility to fit galaxy rotation curves 
TJ" ^ without the need of huge amounts of dark matter [I^, UM, Ej Ej [3 \lM, [3 ■ In addition, it is possible to work out 
■^ , a formal analogy between the corrections to the Newtonian potential and the usually adopted galaxy halo models 
f^ ■ which allow to reproduce dynamics and observations without dark matter [12| . 

OO , However, extending the gravitational Lagrangian could give rise to several problems. These theories could have 

^^ ■ instabilities [I4I , ghost - like behaviors [l^| , and they have to be matched with the low energy limit experiments which 
quite fairly test GR. 

In summary, it seems that the paradigm to adopt /(i?)-gravity leads to interesting results at cosmological, galactic 

and Solar System scales but, up to now, no definite physical criterion has been found to select the final f{R) theory 

2 , (or class of theories) capable of matching the data at all scales. Interesting results have been achieved in this line of 

■ - ' thinking 19, 20, 21, 22] but the approaches are all phenomenological and are not based on some fundamental principle 

as the conservation or the invariance of some quantity or some intrinsic symmetry of the theory. 

In some sense, the situation is similar to that of dark matter: we know very well its effect at large astrophysical 
scales but no final evidence of its existence has been found, up to now, at fundamental level. In the case of f{R)- 
gravity, we know that the paradigm is working: in principle, the missing matter and accelerated cosmic behavior can 
be addressed taking into account gravity (in some extended version), baryons and radiation but we do not know a 
specific criterion to select the final, comprehensive theory. 

In this letter, we want to show that a general exact solution, coming from the request of the existence of a Noether 
symmetry for f{R) cosmological models, matches the two main important requirements that a cosmological solution 
should achieve to agree with data: a transient Friedmann dust-like phase, needed for structure formation, and an 
asymptotic accelerated behavior. Far to be the final model to explain the cosmic speed up, the presence of the Noether 
symmetry could be a physically motivated approach to select viable cosmological models. 
The general features of the theory are the following. Let 

A= [d^xV^fiR)+Ara, (1) 



be the gravitational action where f{R) is a generic function of the Ricci scalar R. GR is recovered in the particular 
case f{R) = —R/16ttG, and Am is the action for a perfect fluid minimally coupled with gravity 
In the metric formalism, this action leads to 4th order differential equations 
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where a subscript R denotes differentiation with respect to R and T™ is the matter ffuid stress-energy tensor. 

In order to derive the cosmological equations in a Friedman-Robertson- Walker (FRW) metric, one can define a 
canonical Lagrangian L ~ L{a,a,R,R), where Q = {a,R\ is the configuration space and TQ = {a,a,R, R} is the 
related tangent bundle on which C is defined. The variable a{t) and R{t) are the scale factor and the Ricci scalar in 
the FRW metric, respectively. One can use the method of the Lagrange multipliers to set i? as a constraint of the 
dynamics. Selecting the suitable Lagrange multiplier and integrating by parts, the Lagrangian C becomes canonical. 
In our case, we have 
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It is straightforward to show that, for f{R) = —R/WttG, one obtains the usual Friedman equations. 

The variation with respect to R of the Lagrange multiplier gives A = /r. Therefore, integrating by parts, the 
point-like FRW Lagrangian is 



C = a^f-fnR) + 6a^fRRRd + 6fRaa^^6kfRa, 



(4) 



which is a canonical function of two coupled fields, R and 
corresponding to the {0, 0}-Einstein equation, is 

Ec ^G fRR. a^ aR + 6 fRaa^ - a^ 



both depending on time t. The total energy Ec, 



if-fRR) + 6kfRa = D. 



(5) 



where D represents the standard amount of dust ffuid as, for example, measured today. The equations of motion for 
a and R are respectively 
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(7) 



where H = a/a is the Hubble parameter. Considering R and a as independent variables, we have, for consistency, 
that R coincides with the definition of the Ricci scalar in the FRW metric. 

The form of the function f{R) and the solution of the system ([5]), (O and ([7]) can be achieved by asking for 
the existence of Noether symmetries. On the other hand, the existence of the Noether symmetries guarantees the 
reduction of dynamics and the eventual solvability of the system [23, [2J, [2q]- Here, we want to seek for viable f{R) 
cosmological models. 

We shall focus our attention on the fact that we need a cosmological solution of the field equations which exhibits 
not only an accelerated phase in recent universe, but also a decelerated period, which lasts for a long time, sufficient 
to allow the formation of structures. This issue has recently been argument of debate since the vali dity of f{R) 
cosmology, which claims to avoid unknown ingredient as dark energy, strictly lies on this possibility [26|. Several 
works on /(i?)-gravity have been devoted to the acceleration and the reconstruction of the models starting from data 
[7|. Numerical treatment is almost obliged and some educated, although arbitrary, guess on the functional form is 
often necessary. On the other hand, /(i?)-cosmology should give rise to standard Friedmann dust-dominated phase, 
which is necessary for the structure formation mechanism, widely accepted and properly working. A first answer 
to this issue was given by means of a numerical reconstruction of the f(R) function [l9|. Here, we want to present 
a general exact solution of the equations, obtained by means of the so called "Noether Symmetry Approach" . A 
summary of the method can be found in [23, l2J, [2^ . 

We ask now for the existence of a vector field 
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such that the Lie derivative of the Lagrangian is zero, i.e. C is conserved and X is a Noether symmetry. It is then 
possible to find 



1/a ; P = -2R/a^ 
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The absolute value is needed, because (with our conventions) we have R < 0. Once the symmetry is found, we 
have an additional constant of the motion, and it is then easy to find a change of variables {a, R} — > {u, u}, such that 
one of the variables is cyclic. We have in fact 



and the new Lagrangian is 



u = a^\R\ ; w = a72 (10) 
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The Noether charge is then Ei ~ u/y/u, leading to immediate integration for u. Introducing the solution into 
Ec — D, and solving for v we obtain 
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\2 



u = -(Sii + Eo) (12) 

288 72 ^48 2j V72Si Ei OEiy " ^ ' 

The parameters Sqj ^i: D, and wq are the integration constants of the equations. They are four since this is a 
general solution of a fourth order problem. 

Coming back to a(i), and setting, for the sake of simplicity a(0) — 0, i.e. vo — 0, we get 



\/a4i'* + ast^ + a2t'^ + ait, (14) 



with 



04 = ; a-i — ; ot = fc ; ai = — — 2A;— — — — -. 
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We see that this solution is a ex t^ for large t, and a oc t^^^, for small t. There is thus room for a smooth transition, 
passing through a period during which the solution approximates reasonably well a Friedmann dust-transient like 
ttf oc t'^'^. In order to see this, we have to consider suitable values of the integration constants a^. All computations 
and results are simplified if we fix the time unit, by setting the current time to = 1. This will not affect the results but 
the value of Hq has to be recast with respect to physical units. We assume also Hq — 1 for simplicity. A unitary value 
for ao can be also set, if no restriction on the value of k is imposed. Finally, we consider a value of the deceleration 
parameter qq = —0.4, which could describe a reasonable current acceleration. These considerations yield a model 
depending only on one parameter. Taking 04 = 0.106, the scale factor turn out to be expressed as 



a= W^[2 + 0.53(t-l)3 + t + 2t2]. (15) 

Comparing this solution with a/ = aoft'^'^ and noting that ao/ must be less than ao, we obtain the very good 
coincidence of Fig[TJ The difference is close to 3% in the interval 2 < z < 4, enough for a phase dominated by 
galaxies. 

It is interesting to come back to the original parameters, in particular for what is concerning the spatial curvature. 
We have k ~ -0.49, which yields ^k^o = kGeff/i^H^al) ~ -0.02, with Ge// = l/[2/'(i?)]. Therefore, this model 
describes a spatially open universe instead of a spatially ffat fc = 0. Indeed, what is physically relevant is not the value 
of k , which is connected with the normalization of a, but the dimensionless parameter fJ^. Moreover, the alleged 
statement fife ~ 0, is obtained from the spectrum of the CMBR radiation and strongly depends on the standard ACDM 
model. Another relevant parameter is the matter content. With our choice of the parameters we get D ~ 0.88, this 
value implies Qm.o — 0.042, which is very close to the expected content of baryonic matter in the Universe. One could 
consider an observer living within a universe described by our model. If this observer is unaware of the fact that the 
fmiction f{R) in the Lagrangian is f{R) — ~ \R\ and not f{R) = R, he would perform all calculations taking into 
account Gn (and not with Geff), obtaining 17^ q ~ 0.29. This value is the expected one for all the matter content in 
the Universe, included the dark matter. Therefore, in this framework, it seems that taking into account dark matter 
could be nothing else but an assumption due to the ignorance of the physical theory behind the cosmological model. 

It can also be noted that ^mfi has nearly the same value of —Q,kfl- Since we have ri,„.o + ^fe.o + ^i?.,o = li the 
current dynamic of this universe results almost totally driven by the curvature, being ^r^ ~ 0.98. 
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Figure 1: Scale factor versus time in standard model (dashed) and our model (continuous). 
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Figure 2: Percentage difference Sa of the two scale factors, for a range in time corresponding to z — 2 -^ 4. It is less than 3%. 



In order to check our model in another way, we consider the distance modulus given by the SNIa and we compare 
our solution with the standard ACDM model, as we know that it fits data very well. Taking as reference the standard 
solution for ACDM model, with flm — 0.27, we get Fig. [3] The coincidence is very good and it is difficult to distinguish 
between the two models. 

Despite these good results, some comments are in order. As we have seen, in our model, the dynamical history of 
the universe is described by the scale factor a{t) ^ t^^^ at early epochs and a{t) ~ t^ at late times giving rise to a 
matter-dust-like stage at intermediate times. This behavior addresses, in principle, the two main issues of dark energy 
models: i) producing a Friedmann-like epoch suitable for LSS formation and ii) an accelerated present epoch stage. 
In Figs, m and [51 we have plotted the behavior of the effective equation of state parameter 
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for our model and compared it with the ACDM model. Clearly, also if the model is accelerating at present epoch 
[z ~ 0), the power is not enough to completely fit the prescription Weff — —1 for the cosmological constant (see 
Fig. [5]). However also the ACDM model does not produce exactly We// = 1 (since there is also the matter component); 
in fact, if we consider Vtm — 0.27, we have We// — —0.73 (in the case ^rn = 0.3, w^ff — —0.7). Therefore, the value 
of our Weff (absolute value) is smaller than the desiderated value, but if we compare this with Weff of the ACDM, it 
is not so far as if we compare it with the pure A-case wa = — 1- 

Furthermore, radiation should be included into dynamics. This fact could destroy the nice feature achieved here, i.e. 
the smooth transition between an unstable dust epoch to a stable, asymptotic accelerated phase. In this perspective, 
more accurate models, including e.g. non-local gravitational corrections, should be taken into account as done in [27| . 




Figure 3: Comparison of the distance modulus S. Our model (continous) and AC DAI (dashed). The agreement is almost 
perfect. 




Figure 4: Weff versus redshift derived from the model. 



Finally, our discussion takes into account only the background while fluctuations are not considered. In fact, at the 
background level, we are able to obtain matter-like regime but things could not work when fluctuations are included 
so one should try to mimic matter - like behavior by modifying gravity or including a dynamical equation of state 
similar to the Chaplygin gas model which well address this goal. This will be the argument of future investigations. 

In summary, we have shown that suitable values of the parameters in the presented general solution (fT4|) allow 
to reproduce the requested behavior of a Friedmann dust - like solution evolving into an accelerated behavior as 
prescribed by observations. This model, physically consistent, has been derived by asking for a Noether symmetry 
in the /(i?) function. The existence of such a symmetry fixes the form of f{R) and allows physically viable models. 
However, starting from this approach, more accurate models should be considered in order to address all the issues 
related to the theory of perturbations and the observational data sets. 
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Figure 5: Comparison of the effective equation of state parameter w^ff for our model (continous) and AC DM (daslied). For 
z larger than 4, the radiation epoch should be carefully considered. 
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